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Abstract 

Differential constraints compatible with the linearized equations of partial differential 
equations are examined. Recursion operators are obtained by integrating the diffe- 
rential constraints. 

One of the standard ways for determining particular solutions to partial differential equa- 
tions is to reduce them to ordinary differential equations which are easier to solve. The 
classical work of Lie about group-invariant solutions generalizes well-known methods for 
finding similarity solutions and other basic reduction methods Bluman and Cole 
proposed a generalization of Lie's method for finding group-invariant solutions, which 
they named the "nonclassical" method. In this approach, one replaces the condition for 
the invariance of the given system of differential equations by the weaker condition for 
the invariance of the combined system consisting of the original differential equations 
along with the equations requiring the group invariance of the solutions. P.J. Olver and 
P. Rosenau proposed a generalization of the nonclassical method |^, They showed 
that many known reduction methods, including the classical and nonclassical methods, 
partial invariance, and separation of variables can be placed into a general framework. 
In their formulation, the original system of partial differential equations can be enlarged 
by appending additional differential constraints (side conditions), such that the resulting 
overdetermined system of partial differential equations satisfy compatibility conditions. 

This work discusses differential constraints compatible with the linearized equations of 
partial differential equations instead of the partial differential equations themselves. The 
relation between differential constraints and recursion operators are examined. For the 
type of equations in the form qt = P{x, t, q, q^, q^x) and qt = P{q, qx, Qxx, Qxxx), recursion 
operators are obtained by integrating the compatible differential constraints. A new type 
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of integrable equations, which are generalizations of the integrable equations of Fokas and 
Svinolupov, are given. Results are also compared with Fokas' generalized symmetry 
and Mikhailov-Shabat-Sokolov's formal symmetry approaches 

We can describe the differential constraint method for evolutionary equations |^] 

qt = P{x,t,q,q^,q^^,...) (1) 

in the following way. First we linearize the given differential equation. In other words, we 
replace q (and its derivatives) in (||) by g + e^' and differentiate both sides of the resulting 
expression with respect to e and take the limit e — > 0, i.e., 

= Dp{^) (2) 

where Dp is the Frechet derivative [Q]. The equation above can also be written as 

where A'" is the order of differential equation, go = Qi = Qx, 92 = Qxx, '^o = ^ , = 
\E'2 = ^xx: and so on. In the classical symmetry approach, (Q) is the main equation, where 
^' is the symmetry of the differential equation, which is a function of x, t, qi. 
The compatible differential constraint is 

= 0, (4) 

where H depends on x, t, qi. If its order (highest derivative in H) is A^, then (^) may be 
written as 

(N-l) 
1=0 

where ^o, Ai, . . . , An-i are functions of x, t, qi. Compatibility of (P) and (§) is given by 
*7V,t - ^'t,7V = 0. (6) 
Using (§) and (^), rewlation (6) leads to 

N-l 

J2 = 0- (7) 

i=0 

Letting 

Wi = 0, (8) 

we obtain a system of partial differential equations among Pi, Ai and their partial deriva- 
tives. The solution of this system will determine the differential constraint (§), which can 
be integrated to give 

= 0, (9) 



where <1> is the recursion operator. 
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Let us now consider differential equations of the following form 
qt = f{x,t,q,q^,q^^). 
The linearized form of (10) can be given as 



(10) 



(11) 



where a, (5, 7 are functions of x, t, q, q^, qxx- We assume a differential constraint having 
the same order as (IC) in the form 



(12) 



where A and B are functions of x, t, q, qx, qxx- The compatibility condition between (|l l\j 

(13) 



and (12) will give the following evolution equations 

At = axx + OLxA + -ixxA + 2-ixAx + ix^ + 27^5 
+2/3^ + Axxl + Axa + 2Ax-iA + 2S^7, 



^4 = laxB + 7^^.S + 27^5^ + -^xAB + 
+ 2Ax-iB + 5^^7 + Bxa. 

The solutions of the system (^)-(|l^) can be given with the linearized equation 



(14) 



= ri'ifxx + 



qx + 2r]r + 7/1 



V {j'qqq^q ^gg) 2 



(Ix + 2r?rgga; 



(15) 



and compatible differential constraint 



ixx Hx' qq 



2 \ 2 
'^qqq'^q ~ '^qqj Qx 



"^^qQx fx 



(16) 



Here ?7i are function of x, t, and r is a function of x, t, q, such that they satisfy 

Ivx-q + rixT = 0, 2rt?7 + ry^r = (17) 



n - ^ (-2r?r7xx + - 27?ir/^ + 4r?r?i,a;) = 0. 



(18) 



It is interesting to note that in this classification the explicit time dependence is crucial 
for obtaining equations (p!7|)-(p!8|), which do not have any derivative with respect to q, 
although r depend explicitly on q. The differential constraint (|l^) can be integrated to 
give 



n I fqqQx , , 7-1 — 1 I r-i — 1 ( 'fqqxfq l^qqfqx 

$ = D + \-r + QxD Vg + qxD ( — ^ 



(19) 



where (Z) ^/) (x) = / /(C) ^C- Because of the last term in (^) and the explicit de- 

J — oo 

pendence of x and t on r, as given by (|T9|), is new recursion operator which is a 
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generalization of the recursion operators obtained by Fokas and Ibragimov-Shabat |1C]. 
The new integrable equation takes the following form 

2 

qt = mxx + + (2r?r + r?i)(?,. (20) 

Although (20) looks like the integrable equations classified by Fokas they are different, 
because of the explicit x, t dependence of r, rj and r/i, which are subject to conditions 
([T7|)-(|l8|). Using transformation q{x,t) — > r{q{x,t), x,t) ^ one obtains 

n = Wxx + (2r?r + ?7i)r^ - ilirq^^q^ + r,^^) - {2r]r + 7/i)r,^ + r,j, (21) 

where r^x, denotes partial derivatives with respect to x, t. However, in the limit 
r{q,x,t) — > r(g), we recover the equations classified by Fokas Ibragimov-Shabat ||l^ 
and Olver |^]. Moreover, equations (^) are more general than the once obtained by 
Svinolupov since he analyzed the case qt = F{x, q, qx, qxx), in which there is no explicit 
time dependence. It is interesting to note that in the limit r{q,x,t) — > r{q,t) (using the 
conditions (17)-(|l8|)), equation (20) will give Svinolupov's equation Q 

qt = r^'qxx + + (2rj'r - + ^) Qcc, (22) 

where r]' and e are functions of t, and r is a function of q and t. 

Equations (0) can be integrated to give r/ = — , where a = a(q). Substitution of r] 
in (|l8|) will give 

ifxxT — ^ar^ + TxT^ir^ + rji^x'f'^ — ftf^ = 0. (23) 

One particular solution can be given by solving the last two terms in the above equation. 
This case {r{q,t)) corresponds to equation (|2^). Equation (|23|) is a coupled second or- 
der partial differential equation and the general solution can be given by the method of 
characteristics [11]. It is hoped that the general solution will appear in a future work. 

Next, we consider differential equations of the following general form, which includes 
the KdV equation, 

qt = P{q,qx,qxx,qxxx)- (24) 
The linearization of (24) takes the form 

'$t = a'^xxx + P^xx + l'^x + 6'^, (25) 

where a, /?, 7, 5 are functions of q, qx, qxx, qxxx- We consider the differential constraint 
having the same order as (|23), i.e., 



^xxx = A^xx + B^x + (26) 

where A, B and C are functions of qx, qxx, qxxx- The compatibility will determine 
algebraic equations 

2n 1 

B = -^, C = --{ax-2aA + 2(3) (27) 
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and evolution equations 

At = a^A + + A^^^rj + 3A^^j]A + 3Alr] + A^a + SA^jjA^, 



1 



= ■^^'^oircxxV - Sa^xV^ - ^o^xAxTj + 2axa 

Pt = ^{^axxxVA - QaxxVA^ - SaxAxxV - l2axAxrjA + Aaxf3 

+PxxxV - Q^xxVA + 4/3^a + 12/3^ ry^2 ^ qa^xPv + 2AAx(3r]A), 

where r] is constant and e is a function of t. 
The hnearized form of the first class is 



Pi ^2 
2 



x+P 



^x + PqQx^, 



with compatible differential constraint 

Q.XX 



XXX — XX 

Qx 



_Pi 2 , ^ 
3r?^^^ 3r/ 



^ _ -"^Pqxx + ^Pqql ^^ 
3r]qx 



where pi is constant and p is a function of q with the condition 
Pm + -^PQ = 0- 



The recursion operator can be obtained by integrating (32). This leads to 



^ = D + — + —Qx - iri^^ + ~^ Pi- 

67] 6rj Srj rj 
The integrable equation is in the form 
Pi s 

qt = mxxx + —q%+pqx- 
6 

The second class is given by the linearized equation 
= V^xxx + [eiql + 2eie2qx + 263] ^ 
with compatible differential constraint 

qxx 



X + £2 



£]_ 2 _^ 2eie2 ^ 2e3 
3?7 ^ 3r] ^ 3r] 



(e^ei - 263) 
3r]{qx + €2) 



(28) 
(29) 

(30) 



where ei, €2, £3 are constants. The recursion operator can be obtained by integrating 
which gives 



-.2 2e3 €2 2 , 2eie2 

1 c 

3?7 3r/ 



<^> = D' + ^ + + =^qx -^iqx + e2)D-\x 

6rj 6r] 



The integrable equation is in the form 

ei 3 €162 2 , 
qt = mxxx + -^qt + —1x + ^sqx- 



(31) 
(32) 
(33) 

(34) 

(35) 
(36) 

(37) 
(©' 

(38) 
(39) 
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Equations (|35|) , (l39|) were classified by Fokas [Q, Ibragimov-Shabat |10|]. We obtained the 
recursion operators (|3^), ( |38| ) by integrating the differential constraints. 
The third type of equation is given by the linearized equation 



3 3 



(40) 



and differential constraint 

nxxx + Aife + X2qx 



+ 



Qxx + AiQx + hq 

^3{Qxxx + ^iQxx + ^2Qx) 



qxx + Aig^ + A2g 
The recursion operator is 
$ = Z)2 + AaZ) + A4. 
The integrable equation is given as 

qt 



^4:{qxxx + ^iqxx + A2fe) 

+ ^iqx + A2q' 



(41) 



Asg + ^^sAsfe + -^^A^^qx + Aeg, 



(42) 



(43) 



where Ai, A2, A3, A4, A5 and Ag are constants. Note that Rabelo and Tanenblat also 
obtained linear equation using the classification method of pseudo-spherical surfaces with 
Gaussian curvature (—1) []T2| . 

According to Fokas, an integrable equation has infinitely many generalized symmetries 
which are the solutions of (§) . The existence of generalized symmetry manifest itself by the 
existence of an admissible Lie-Backlund operator. The existence of infinitely many sym- 
metries is expressed by the existence of a recursion operator. There is a close relationship 
between Lie-Backlund operator and linearized equation H]. Because, if Fokas' admissible 
Lie-Backlund operator is applied on the evolution equation, we obtain the Frechet deriva- 
tive of our differential equation under consideration or linearized form of our differential 
equation. Recursion operators, in our method, are obtained by the integration of the 
compatible differential constraints. 

Let us briefly recall Olver-Fokas symmetry approach and the differential constraint 
test: 

Olver-Fokas symmetry test: The equation qt = P[q] is integrable if there exists 
infinitely many non-Lie point symmetries, or equivalently, one non-Lie point sym- 
metry and a recursion operator. The recursion operator and the time-independent 
part of the linearized differential equation form a Lax pair + Dp] = 0. 

Differential constraint test: The equation qt = P[q] is integrable if there exists 
a differential constraint = 0, such that it is compatible with the linearized 
equation = Dp{^)- The compatibility condition is Ht + [H,Dp] = 0. 

There are several methods to examine the integrability of nonlinear partial differential 



equations [13|, although in two dimensions most of these methods imply each other. 
Infinite sets of conservation laws, infinite number of symmetries, and the bi-Hamiltonian 
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structure, are a number of remarkable properties, to name but three. The recursion oper- 
ator plays an important role in the formulation of these recursive properties. Firstly, the 
family of integrable equations can be written in a compact form using recursion operators. 
The other role of recursion operator is associated with the Hamiltonian treatment of inte- 
grable equations. Recursion operators determine Hamiltonian structures through certain 
Poisson brackets [0]. Further analysis of the recursion operators leads to the concept of 
bi-Hamiltonian structures, which can be given by factorizing the recursion operator. 

The main idea in this work is to give a new definition of integrability. A partial 
differential equation is integrable if its linearized equation is compatible with a differential 
constraint. Using differential constraints compatible with the linearized equations, we also 
give an answer to the question "The deep connection between the direct reduction and 
recursion operators" of Olver [^. 
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